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Abstract 

We show that the hmiting eigenvalue distribution of random symmetric ToepUtz 
matrices is absolutely continuous w it h den sity bounded by 8, partially answering a 
question of iBrvc. Dembo and Jiang The main tool used in the proof is a 



spectral averaging technique from the theory of random Schrodinger operators. The 
similar question for Hankel matrices remains open. 

1 Introduction 

An n X n symmetric random Toeplitz matrix is given by 

T„ = ((ct|j-fc|))o<i,fc<n 

where (oj)j>o is a sequence of i.i.d. random variables with Var(ao) = 1. For a m x m 
Hermitian matrix A, we denote by 



^ m 



m 



the empirical eigenval ue distribution of A, where A-, -, 1 < j < m are the eigenvalues of A 



counting multiplicity. iBryc. Dembo and Jiang! (120061 ) established using method of moments 



that with probability 1, /i(n~^/^T„) converges weakly as n — t- oo to a nonrandom symmetric 
probability measure 7 which does not depend on the distribution of Oq, and has unbounded 
support. They conjecture (see Remark 1.1 there) that 7 has a smooth density. In this note, 
we give a partial solution: 

Theorem 1. The measure 7 is absolutely continuous with density bounded by 8. 
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The actual bound we get is = 7.20 . . ., but we do not expect it to be optimal. 

It seems that the method of moments is of little use in determining the existence of 
the absolute continuity of the limiting eigenvalue distribution. Indeed our proof goes along 
a completely different path. We make use of the fact that the spectrum of the Gaussian 
Toeplitz matrix can be realized as that of some diagonal matrix consisting of independent 
Gaussians coniugated by an appropriate projection matrix - a fact observed in a recent paper 



Sen and Virag] (120111 ) . The next key in gredient of our proof is a spec t ral av eraging technique 



(Wegner type estimate) developed by I Combes. Hislop and Mourrd (119961 ) in connection to 
the problem of localization for certain families of random Schrodinger operators. 

Our proof does not establish further smoothness property of 7. The absolute continuity 
for the limiting distribution of random Hankel matrices also remains open. 



2 Connection between Toeplitz and circulant matrices 

Since 7 does not depend on the distribution of Oq, from now on, we will assume, without 
any loss, that (aj)j>o are i.i.d. standard Gaussian random variables. The remainder of the 
section we rec all some facts about t he connection between Toeplitz matrices and circulant 



matrices from ISen and Viragi (I2OIII ). Let T° be the symmetric Toeplitz matrix which has 
V^oo on its di agonal i nstead of oq. It can be easily shown (e.g. using Hoffman- Wielandt 
inequality, see bhatial ^99^) that this modification has no effect as far as the limiting 
eigenvalue distribution is concerned. 

T° is the nxn principal submatrix of a 2nx2n circulant matrix C2n = mod 2n)o<i,j<2n-i, 
where bj = aj for < j < n and bj = a2n-j for n < j < 2n, bo = -\/2ao, bn = V^an- In other 
words, 

where Q2n = 

0„ 0„j \0n Or. 

The circulant matrix can be easily diagonalized as (2n)^^/'^C2n = U2nD2„U2„ where U2n is 
the discrete Fourier transform, i.e. a unitary matrix given by 

/ 2iTijk\ 



Q2nC2nQ 



2n 



(1) 



U2n(j, k) 



exp 



[ 2n ) 



,0 <j,k <2n-l 



and D 



2n 



diag(c?o, rfi, . . . , rf2n-i), where 

2n-l 



k=0 



feexp 



27rij k 
2n 



n-l 



V2ao + V2an + 2 ^ cos 



k=l 



f27rjk\ 
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Clearly, dj = d2n-j for all n < j < 2n. Also, ('ij)o<i<n are independent mean zero Gaussian 
random variables with Var((ij) = 1 if < j < n and Var((i_,) = 2 if j G {0,n}. Define 

P2n := U;nQ2nU2n SO that 



Check that is a Hermitian projection matrix with P2n(j, j) = 1/2 for all j. For notational 
simplification, we will drop the subscript 2n from the relevant matrices unless we want to 
emphasize the dependence on n. 

3 Proof of the main theorem 

For a vector u G C™, let cr(A, u) be the spectral measure of matrix A at u. For a finite 
measure u on M, its Cauchy-Stieltjes transform is given by 



Let ¥,fi(n ^^^T° ) denote the expected empirical eigenvalue distribution of n ^^^T° which is 
defined by En{n-^/'^T°J{B) = E[fx{n-^/^T°^){B)] for all Borel sets B. 

Lemma 2. Let (ej)o<j<2n-i be the coordinate vectors o/R^". Then 

2n-l 

^ i=o 

Before we start proving the above lemma, we state a simple fact about spectral measures 
of Hermitian matrices without proof. 

Lemma 3. Let A be an m x m Hermitian matrix. Let Ui, U2, . . . , and Vi, V2, . . . , be 
vectors in satisfying Yli=i ^i^i = Sj=i ^i^j ■ Then 



{2n) "'^^^U2„Q2nC2nQ2nU2n — P2nD2nP2n- 



(2) 




k e. 



^a(A,uO = ^a(A,v,). 



1=1 j=i 



Proof of Lemma [H By ([T]) , we have 




j=0 
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Changing basis as in ([2]), we can rewrite this as 

Pi '^^^ Pi "'^^ 

— V(U*e„ (PDP - zI)-^\J*e.) = — V s(z; (xfPDP, U*e,-)). 
n ^-^ n ^-^ 

j=0 j=0 

Now by Lemma |3] and the fact that ^"=0 U*eje*U = X]j=o^ P^j^jP) deduce 
s{z- f^in-'/'TD) = ^ ^(Pe,, (PDP - ziy^Pe,). 



(3) 



j=Q 



The lemma now follows by taking expectation on both sides of ([3]) and by observing that 
for a fixed z G C, Im(2;) 7^ 0, the map u s{z] v) is linear and hence commutes with the 
expectation. □ 

Next we will prove a key lemma about the uniform bound on the Stieltjes transform of 
the expected empirical eigenvalue distribution of Toeplitz matrices. 

Lemma 4. For all n, we have 

sup |s(z;E/i(n-^/2T°))| < 16^2. 

z:Im(2;)>0 

The above lemma will be a direct consequence of the following result of 



Combes et al. 



( 19961 ) on the spectral averaging for one parameter family self-adjoining operators. 



Proposition 5 (ICombes et al.l ( 1l996l )). Let Hx,X E M be a C'^-family of self-adjoint op- 
erators such that D{H\) = -Do C "H VA G M, and such that {H\ — z)~^ is twice strongly 
differentiahle in A for all z, lm.{z) 7^ 0. Assume that there exist a finite positive constant Cq, 
and a positive hounded self-adjoint operator B such that, on Dq 



Also assume H\ is linear in X, i.e., H\ 



dX 



>coB' 



dA2 



(4) 



0. Then for all E E M and twice 



continuously differentiahle function g such that g,g',g" G //""^(M) and for all (f eH, 



sup 

<5>0 



/ g{X){^,B{Hx-E-t6)-'Bip)dX 
Jr 



<Co'(ll^7l|i 



1^ 111 



\9 \\i)\m 



(5) 



Proposition [5] is an immediate corollary of Theorem 1.1 of ICombes et al.l (119961 ) where 
instead of Hx = 0, it was assumed that \Hx\ < CiHx- The vanishing second derivative 
assumption shortens the the proof by a considerable amount. We have included a proof of 
the above proposition in the appendix to make this paper self-contained and also to make 
constant in the bound ([5]) explicit. 
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Proof of Lemma\^ Set Ej = e^e* + e2„-je2„_^- for I < j < n, and = e^e* for j G {0, n}. 
Take 

Bj = Peje*P or Pe2n-je*2„_jP for 1 < j < n and = Pe^e^P for j e {0, n}. (6) 

Fix < j < n. We apply Theorem [5] with Hx = P(D + (A — (ij)Ej)P. In words, we 
replace dj and d2n-j by A in PDP to get H\. Note that Hx is random self-adjoint operator 
which is a function of {dk : < k < n,k ^ j}. Also, Hx is linear in A and so, Hx = 0. 
Since Hx = PE^P > B^- = P(j,j)"^B^, the condition is satisfied with B = Bj and 
Co = 2 since P{j,j) = 1/2. Take g = (pj where (pj be the density of Z for < j < n 
or the density of \/2Z for j G {0,^}, Z being a standard Gaussian random variable. It 
is easy to check that \\g\\i = 1, \\g'\\i < \\9"\\i ^ 2. Then plugging (p = e.j or e2„-j 
and Bj = Peje*P or Pe2n~j^2n-j^ and taking expectation w.r.t. the remaining 
randomness {dk : < k < n, k ^ j} , we obtain 

sup P(j,jf E(Pe,-,(PDP-zI)-iPe,-) < Cq i(||(7||i + + ||/||i) < 2. (7) 

z:lm(z)>0 

The lemma is now immediate from and Lemma [2j □ 

Proof of TheoremUl It follows from the inversion formula, z/{(x, y)} = lim^j^o " /J Im(s(i? + 
i6; v))dE for all x < ?/ continuity points of v that if for some probability measure /i, 
sup^.jjjj(2)>o Ini(s(z; < K then is absolutely continuous w.r.t. Lebesgue measure and 
its density is bounded by -k^^K. 

Note that 5(2;; Eyu(n^-'^/^T° )) — i- s{z \ 7) as n — > ■ co f or each z G C, Im(2;) > since 



E/i(n ^^^T°) converges weakly to 7 (see lBryc et al.l (I2OO6I )). So by Lemma HI it follows that 



sup \s{z\ 7)1 < 16^2 < Stt 

^:Im(2)>0 

which completes the proof of the theorem. □ 



Appendix 

Proof of Proposition Define for e > and < 5 < 1 , 

R{\e,5):={Hx- E + i5 + itHx)-^ (8) 

and set 

K{X,e,6) := BR{X,e,6)B. (9) 
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Note that from assumption (jlj) , 

-Im(^, K{X, e, 6)^) = BR{X, e, + eHx)RiX, e, 6)B^) > Coe\\K{X, e, 6M\ 
which, coupled with Cauchy-Schwarz inequahty, imphes that G "H, \\(p\\ = 1, 
\\K{X, e, 5M > -Im(^, K{X, e, 5)ip) > coe||K(A, e, SM'. 

Now define 

F{e,5) := / g{X){^,KiX,e,S)^)dX. 
Jr 

Inequahty (fTU]) imphes the bound 

F{e,6)<{eco)-%\\i. 

Now differentiating F w.r.t. e, we obtain 
.(iF(e,5) 



(10) 



de 



g{X) {if, BR{X, e, 6)HxRiX, e, 5)B^)dX 



d 



- / ^/(A) — (^,ir(A,e,5)^)dA. 



where the last equality follows from the fact Hx = 0. Therefore, from ffTOj) and by integration 
of parts, 

dF{e,6) 



de 



g'{X){^,K{X,e,S)y^)dX 



< (eco)-'lklli. 



(12) 



By integrating the differential inequality ( IT2|) and using the bound (fTTj) . we can improve the 
bound for F as 

|F(e,5)| < coi^^'lli ■ |loge| + \F{1,5)\ < c^%'h ■ | loge| + Coi^^Hi, Ve G (0,1). (13) 

Now if we consider the function F{e,6) := f^g'{X){ip, K{X,e,6)ip)dX, then by replacing the 
function g by its derivative g' in (fT5]) . we deduce that 

|^(e, 5)1 < Co i/lli ■ I logel + c,%'h, Ve G (0, 1) 

which further implies that 
dF{e,6) 



de 



<coi/||i-|loge|+co-l(7'||i, VeG(0,l). (14) 
Again integrating dH]), we get 

\F{e,S)\ < c,\\\g"h + WWi) + \F{l,S)\ < c,\\\g"\U + H^^'IK + ll^?lli), (15) 

which holds for all e, 5 G (0, 1). The proof of the Proposition now follows from the fact that 
i?(A, e, 6) converges weakly to {Hx — E + i6)~^ as e ^ 0+ provided 6 > 0, and the dominated 
convergence theorem since \J^g{X){(p, K{X, e, 5)(p)dX\ < C, by f|T5l) . □ 
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